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Abbildung 5: Material- und Energietransport in WECHSL 
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Abbildunq 8: 


Schmelzenkaverne im Reaktorfundament nach 42 h 
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Einleitung 





Bei der Analyse von Kühlmittelveriuststörfällen von Leichtwasser- 
reaktoren ist der Massenstrom des Kühlmittels aus einer Bruchstelle 
Das Kühlmittel strómt dabei sich zeitlich 


stark andernd als zweiphasiges Gemisch aus Wasser und Wasserdampf 
aus der Bruchstelle. 


eine maßgebliche Größe. 





Weltweit werden derzeit technische Großversuche zur Simulation 
von Kühlmittelverluststörfällen durchgeführt; z.B. 
HDR-Programmes - BRD, 


im Rahmen des 
bei den Flutexperimenten in der Large Scale 
Reflood Test Facility - Japan, bei den Blowdown Untersuchuncen 
LOBI - Italien und im LOFT-Programm - USA. 





Wenn die Übertragung der Ergebnisse aus diesen Simulationsver- 
suchen auf reaktortypische Verhältnisse gelingen soll, 
die Zweiphasenmassenströme mit hoher Genauigkeit ermittelt 


dann müssen 




























Die Messung 


von zweiphasigen Massenstrómen ist gegenüber der 
einphasigen Massenstrómen dadurch erschwert, daß 


und Geschwindigkeitsverteilung im Rohrquerschnitt 
Formen besitzen kónnen. 


Messung von 





die Phasen- 
vielfältige 
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Die Vielfal: 
Weise durch 


der Massenverteilungen wird in vereinfachender 
typische Strömungsformen charakterisiert. 
zeigt die bei waagerechter Rohrstrómung charakteristischen 





Strómungsformen: Disperse Blasenstrómung, Blasenstrómung, 
Schwallstrómung, Ring-Spritzerstrómung, Schichtenstrómuna und 
Wellenstrómung. Welche Strómungsform sich einstellt ist iw 
wesentlichen abhängig von den Massenströmen der einzelnen Phasen 
und den Stoffeigenschaften. Man versucht die Grenzen zwischen den 
einzelnen Strómungsformen in sog. Strómungskarten darzustelien. 
In Abb.1 ist als Beispiel dafür die Strómungskarte von Mandhane 2 
al. /1/ enthalten, die als Koordinaten die bezogene Gas- und 
Flüssigkeitsgeschwindigkeit (superficial velocities) verwendet. 


Abb. 1 macht deutlich, da8 solche komplexe Strómungen hohe An- 
forderungen an ein universell einsetzbares Zweiphasenmassenstrom- 
Me8gerát stellen. Die Notwendigkeit, die Zweiphasenmassenstróme 
möglichst genau messen zu können, hat deshalb in den vergangenen 
Jahren bei verschiedenen Forschungsinstitutionen zur Entwicklunc 
untersch:edlicher Meßverfahren geführt. Ein Überblick über der- 
artige MeBverfahren ist in /2-4/ enthalten. 


Auf Anregung des Sachverständigenkreises Notkühlung im Jahre 1974 


wurden im Rahmen des Projektes Nukleare Sicherheit im Institut für 


Reaktorbauelemente Versuchseinrichtungen für stationáre Luft- 
wasser- sowie Dampf-Wasser-Strómung geschaffen. Mit Hilfe dieser 
Versuchskreisiäufe, die im Detail in /5/ beschrieben sind, werden 
im IRB seit einigen Jahren Arbeiten durchgeführt, die als Ziel 
haben, verschiedene Zweiphasenmassenstrom-Messverfahren unter 
vergleichharen Bedingungen zu erproben, zu kalibrieren und die 
Auswerteverfahren methodisch zu verbessern. 


Im Rahmen dieser Zielsetzung besteht eine enge Kooperation mit 
Forschungsqruppen anderer lIastitutionen, die Meßverfahren meist 
für spezielle Sicherheitsexperimente entwickeln. Tabelie 1 zeint 
eine Zusamaenstellung der im IRB getesteten Instrumentierungen 
in der zeitlichen Reihenfolge des Einsatzes. 


Im folgenden werden einige dieser Verfahren nerausoegriffen und 
näher erläutert. Dabei wird die Darstellun4 der MeBverfahren 
nach einer gewissen Rangfolge (Klassifizierung) der Verfahren 


vorgenommen. Eine Beyründung für diese Rangfolge folgt später. 





2. Beschreibung der Meßverfahren 





Das Radionuklidverfahren />/, /7/ basiert im wesentlichen auf 





einer Laufzeitmessung von Strómungsanzeigern und einer Messung 

des Gasvolumenanteils mit Hilfe eines Gamma-Densitometers (Abb.2). 
Die Strómungsanzeiger bestehen aus kleinen Mengen von kurzlebigen 
Radioisotopen. Die flüssige und die gasfórmige Phase werden dabei 
getrennt markiert. In die Fliissigkeit wird das Radioisotop Mn 56 
als Mangan-Azetatlósung eirgebracht. Das Gas wird duch das \rgon- 
Isotop 41 markiert. Eine hohe Verweilwahrscheinlichkeit der Iso- 
topensubstanzen ist nur in den jeweiligen Phasen vorhanden.Stromab der 
Injektionsstelle werden die von der Zweiphasenströmung transportier- 
ten Isotopenwolken durch Szintillationszähler in zwei Meßebenen 
detektiert. Weil die Gamma-Energien der beiden Isotope unterschied- 
lich groß sind. können die Isotope mit elektronischen Energiefenstern 
getrennt detektiert werden. Eine Korrelation zwischen den Detekti- 
onsimpulsen ergibt eine Laufzeit. Aus dem räumlichen Abstand 

zwischen der Positionen der Detektoren können die Phasengeschwindig- 
keiten des Gases Teo, und der Flüssiqkeit Ya] errechnet werden. 

(Der zusätzliche Index R soll darauf hinweisen, daß die Phasenge- 
schwindigkeiten mit der Radiotracer-Technik ermittelt wurden). Die 
Massenstromdichte G (Verháltnis von Massenstrom ñ zu Rohrquer- 
schnitt A) wird mit folgender Beziehung ermittelt 


Rg (i-a Ip, V 


Rl 


wobei % und ገ die Dichten der Gas- und Fijssiqphase sind und 
der Gasvolumenanteil (void fraction) ist. Der Gasvolumenantei | 
wird üblicherweise mit Hilfe eines Gamma-Densitometers gemessen. 
Dieses Instrument basiert auf der Schwächung der Intensität eines 
dünnen Gammastrahls beim Durchdringen von Materie. Die Messung 
ergibt eine entlang des Strahls gemitteite Dichte. Um den über den 
Rohrquerschnitt gemittelten Wert zu bestimmen, werden in der 
Praxis mehrere Gammastrahlen verwendet. Zwischen der von einem 
aamma-Densitometer gemessenen mittleren Dichte c, und dem mittle- 


ren Damofvolumenanteil besteht folgender Zusammenhang 


rür die in dieser Arbeit behandelten Ergebnisse wurde das 3-Strahl 

Gamma-Densitometer verwendet, das ein Teil der spáter beschriebener 
LOFT-Instrumentierung ist. In der Zwischenzeit wurde im Rahmen der 

Entwicklung des Radionuklid-Verfahrens ein Gamma-Densitometer ent- 

wicxelt, das in /7/ und /8/ beschrieben ist. 


Das Rad.onuklid-Verfahren wurde bisher bei den Marviken I und II- 
Experimenten (Schweden) und einigen HOR-Experimenten (DPRD) eingesetzt. 


Ein Einsatz in den LOBI Experimenten (Euratom Ispra, Italien) wird 
derzeit vorbereitet. 


Das True Mass Flow Meter (TMFM) /9/, /10/ entspricht im Aufbau einer 
radialen Strómungsmaschine (Abb. 3): Das axial strómende Zweiphasen- 
gemisch wird durch einen mit konstanter Drehzahl umlaufenden Rotor in 
radiale Richtung umgelenkt. Die umgelenkte drallbehaftete Strömung 

übt dann auf ein Leitschaufelsystem ein Moment aus. Dieses Moment wird 
mit DehnungsmeBstreifen gemessen und ist ein direktes Maß für den 
gesamten Massenstrom bzw. die gesamte Massenstromdichte. Der 
Zusammenhang zwischen Massenstromdichte G und Moment M ist durch 
folgende Gleichung gegeben: 
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"men j Au r* 

mit w= Winkelgeschwindigkeit, r = Schaufelradius. In dieser Beziehunq 
ist allerdings die Annahme enthalten, dai der Schaufelwirkunasarad 
identisch 1 ist. Dies trifft nur für eine ideale Strómungsmaschine 
zu. Die Größe des t2tsächlich vorhandenen Schaufelwirkungsarades 

kann mit empirischen Beziehungen aus dem Strómungsmaschinenbau 


abgeschátzt oder durch Messungen in reiner dasserstrómung ermittelt 
werden. 


Das in dieser Arbeit betrachtete TMFi! wurde für das PAS-Cosima-Vor- 
haben entwickeit. Ergebnisse, die im Rahmen dieses Vorhabens in in- 


stationärer Strömung erzielt wurden, sind in /11/ und /12/ entsal- 
ten 


_ GE — 


Die MeBinstrumentierung, die in Abb. 4 dargestellt ist, besteht 
aus einer Stauscheibe (Drag Disc), einem Turbinenradzähler 
(Turbine Flow Meter) und einem 3-Strahl-Gamma Densitometer. Diese 
Art der Instrumentierung wird in mehreren Sicherheitsexperimenten 
verwendet, wobei die konstruktive Gestaltung der Einzelinstrumente 
verschieden ist. In diesem Beitrag wird die Instrumentierung 
betrachtet, die für LOFT (Loss of Fluid Test) in Idaho, USA /13/ 
verwendet wird. Hier ist die Stauscheibe und der Turbinenradzáhler 
in einer Einheit zusammengefaßt, dem sogenannten Drag Disc Turbine 
Transducer (DTT) /14/. Der Durchmesser des DTT ist kleiner als 

der Rohrdurchmesser; die Messung erfolgt daher in gewissem Maße 
lokal (free field). Es besteht die Frage, auf welche Strömungs- 
eigenschaft die einzelnen Meßinstrumente ansprechen. Es wird 
angenommen, daß die Stauscheibe die querschnittsgemittelte Im- 
pulskraft mißt, d.h. die folgende Beziehung gilt: 
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2 B s 
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Für den Turbinenradzähler wird angenommen, daB das Ausgangs- 
signal ein Ma3 für den gesamten Volumenstrom ist. Dividiert man 
den Volumenstrom durch díe Rohrquerschnittsfláche, so erhält man 
eine Geschwindigkeit, die gleich der homogenen Geschwindigkeit 


V nom 151: 


w, = av +(1-a)V, (5) 


Das Gamma Densitometer 511] die mittlere Dichte des Gemisches 
messen, die folgendermaßen gebildet ist: 


o,” ap, *(1-0)0, (6) 


Aufgrund der 3 Strahlen, die über den gesamten Rohrquerschnitt 
verteilt sind, wird mit dem Gamma Densitometer in sehr víel 
höherem Made ein querschnittsgemittelter Wert bestimmt, als 
mit den DTT-Signalen. 
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Die Berechnung der Massenstromdichte erfolgt oft nach den fol- 
genden Beziehungen, wobei jeweils die Meßwerte von zwei MeBgeršten 
kombiniert werden: 


(7) 
(8) 


(9) 


Jiese Auswertebeziehungen sind zwar sehr einfach, sie gelten 
jedoch, wie sich leicht zeigen läßt (vgl. z.B. /14/), nur für 
homogene Strómung. 


Eine andere Möglichkeit, die Massenstromdichte zu bestimmen, 
besteht darin, mit den 3 unabhängigen Meßwerten (ov^) pg» ^, 

und Y, mit Hilfe von geeigneten Modellen den Dampfvolumenantei | 
a und die Phasengeschwindigkeiten Ya und Y, zu ermitteln (vgl. 
z.B. /15/). Für die Massenstromdichte wird dann die Beziehung 
verwendet, die ohne Einschrinkung gültig ist 


G = ap V * (i72) 9. Vi (10) 


Eine weitere Meßmethode, die zum Vergleich herangezogen wird, 
besteht aus einer Verturidüse, einem Turbinenradzähler und einem 

1 Strahl Gamma-Densitometer. (Abb. 5). Im Gegensatz zum LOFT-MeBein 
satz überstreicht hier der Turbinenradzähler den gesamten Rohrquer- 
schnitt (full flow). Dieser Meßeinsatz wurde in CEN Grenoble entwickelt 
/16/ und wird bei den PHEBUS Experimenten in Cadarache, Frankreich 
eingesetzt. In diesem Beitrag wird zur Bestimmung der Massenstrom- 
dichte nur die Druckdifferenz über die Venturidiise ^p und die 
Drehfrequenz des Turbinenradzühlers f verwendet: 

















Gy - Ky Í oy(ap+gh(o,-0,)) 
Gr = Kr ፍገ 


wobei Ky - Beiwerte sind, die in einphasiger Strömung best: wet 
werden, g h(0,70,]) ist ein Korrekturterm, der bei vertikalem 
Einbau der Düse auftritt. 


lur Bestiamung der Massenstromdichte aus den obigen Beziehungen 
muB im wesentlichen die Annahme gemacht werden, daß Be = Pre 
eine Annahme, die a priori nicht erfüllt ist. 
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Gülti 
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Ein Bei 
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dichte 
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(z.B. ! 
sitome 
Strómu 


ilung der MeBverfahren 





2 werden im folgenden analogy zu der Vorgehensweise 
ie Meßsysteme nach der Art der Massenstrombestimmung ein- 


und für jede Gruppe ein Beispiel angeführt: 


ahren wird als vollständig bezeichnet, wenn 
‚swertebeziehung ohne Einschränkung gültig ist und 
der Auswertebeziehung vorhandenen Größen direkt 


sen werden. 


dingungen sind für das TMFM und das Radionuklid-Ver- 
erfüllt. 


verfahren wird als unvollständiges Verfahren bezeichnet, 


tweder 
uswertebeziehung nur für spezielle Strómungsformen 


gkeit hat 


eßwerte über Modelle mit den eigentlichen GrüBen ver- 
t sind, wobei der Gültigkeitsbereich der Modelle 


nicht genau bekannt ist. 


spiel fur diese Gruppe ist die Massenstrombestimmung 

| 3 MeBwerten von Stauscheibe, Turbinenradzähler und 
lensitometer. Hier ist die Sezienung für die M:ssenstrom- 
zwar ohne Einschränkung gültig; die Modellbeziehungen 

| Stauscheibe und den Turtinenradzáhler sind jedoch un- 
jigend, wie die praktische Erfahrung zeigt (vgl. /15/). 


Gruppe der unvollständigen Verfahren gehürt ebenfalls 
ssenstrombestimmung mit 2 der oben genannten 3 MeBwerte 
sie in Tab. 2 dargestellt, díe Kombination aus Gamma-Den- 


ter und Stauscheibe). Hier wird implizit eine homogene 


ng angenommen. 














Entsprechendes gilt für die Massenstrombestimmung mit den Meß- 
werten der Venturidüse und des Turbinenradzählers: Hier wird 

die Dichte nicht direkt gemessen, sondern bei der Ermittlung 

der Massenstromdichte eliminiert, was nur zulássio 1st, wenn 

die "Modellvorstellung" richtig ist, da3 die für beide Instrumente 
maßgebenden Dichten gleich sind. 


Als Konsequenz der in Tab. 2 getroffenen Einteilung der Meß- 
verfahren kann man sagen, daB bei den vollständigen Verfahren 

auf eine Kalibrierung in zweiphasiger Strómung im wesentlichen 
verzichtet werden kann, die MeBbereichsgrenzen und Genauigkeiten 
können aus Messungen in einphasiger Strömung in befriedigender 
Weise bestimmt werden. Bei den unvollständigen Verfahren ist eine 
Kalibrierung in zweiphasiger Strömung unbedingt erforderlich. 


4. Dampf-Wasser-Kreislauf 





Eine ausführliche Darstellung der Versuchskreisläufe ist in /5/ 
enthalten, hier wird nur kurz der Dampf-Wasser-Kreislauf beschrie- 
ben (Abb.6): Zwei Dampferzeuger, von denen einer leicht über- 
hitzten Dampf, der andere leicht unterkühltes Masser bereitsteilt, 
versorgen den Versuchsstand. Stromab der Mischkammer befindet 

sich die Teststrecke, die mit den zu untersuchenden Instrumen- 


ten bestückt ist. 


Die Teststrecke ist mit einem 5-Strahl-Gamma-Densitometer /17/ 
(Teststreckendurchmesser 50 mm) und mit lokal messenden Impedar:- 
Sonden /18/, /19/ ausgestattet. Diese MeBinstrumente sind Teil der 
allgemeinen Kreislaufinstrumentierung und haben den Zweck, die 
Strömungszustände der Zweip.asenströmung in der Teststrecke zu 
identifizieren /20/, /21/. 


In einem Kondensator wird das Zweiphasengemisch kondensiert und 
über Pumpen in die Dampferzeuger zuriickgefirdert. 







Wesentliche Bestaniteile des Kreislaufes sind Meßblenden, 
welche jeweils getrennt den einphasigen Massenstrom des Wasser- 
dampfes und des Wassers vcr dem Eintritt in die Mischkammer er- 











mitteln. Die Summe der Massenstróme ergibt den gesanten Massen- 
strom, der als Referenzwert zum Vergleich mit den Werten benutzt 
wird, die von den eingesetzten MeBverfahren gemessen werden. 


5. Ergebnisse 


Im folgenden werden die Mefergebnisse der einz^lnen Verfahren 
miteinander verglichen. Die Tests wurden in Dampf -Hasser-Strómung 
durchgefiihrt bei vergleichbarer Variation von Massenstrom, Dampf- 
gehalt und Druck. Die Versuchsmatrix umfaßte folgende Strómungs- 
formen: Weilenstrómung, Schwallstrómung und exzentrische Ring- 
Spritzerstrómung, also Strümungsformen, bei denen ausgeprágte 
Phasen-und Geschwindigkeitsunterschiede im Rohrquerschnitt auf- 
treten. Das Radionuklid-Verfahren und der LOFT-MeBeinsatz wurden 
gemeinsam in waagerechter Strümung getestet. Das TMFM, dessen 
Rohrachse senkrecht orientiert sein muB, befand sich am Ende einer 
horizontalen Teststrecke. Der PHEBUS-MeBeinsatz wurde zusammen mit 
der bei den PHEBUS Experimenten vorhandenen Rohranordnung (vgl. 


Abb.5) getestet. 


In den nächsten Abbildungen ist jeweils die von einem getesteten 
Verfahren gemessene Massenstromdichte über der Referenzmassen- 


stromdichte aufgetragen. 


Abb. 7 zeigt die TMFM-Ergebnisse: Die Streuung der Versuchs- 
punkte ist sehr klein, der mittlere Fehler des Verhältnisses 
TMFM-Massenstromdichte zu Referenzmassenstromdichte beträgt 

= 3 x, die Standardabweichung =10 í. Die Abweichungen von den 
Referenzwerten sind praktisch unabhängig von der Massenstrom- 
dichte und, was aus der Abb. 7 nicht zu ersehen ist, von Dampf- 
gehalt und Druck, d.h. praktisch unabhängig von der Strömungs- 
form in der Teststrecke /22/. vie Abweichungen sind im wesent- 
lichen dadurch bedingt, daß in der Auswertung ein Schaufelwir- 


kungsgrad von 1 angenommen wurde. 


Abb. 8 zeigt die Ergebnisse, gewonnen mit dem Radionuklid-Ver- 
fahren (/14, /23/): 


s/ By ‘wosjsuassey 


Auch hier ist eine gute Obereinstimmung zwischen gemessenen 
und Referenzwerten vorhanden. Der mittlere Fehler der Verhält- 
nisse der Massenstromdichten beträgt + 5 5, die Standardabwei- 
chung 7 %. Die größten absolsten Abweichungen treten im 

oberen Bereich der Massenstromdichte auf. Es wurde in 

/15/ gezeigt, daB sie im wesentlichen durch die Bestimmung 

des querschnittsgemittelten Dampfvolumenanteils aus den 
MeBwerten der einzelnen Gamma-Strahlen verursacht werden. 
Verbesserungen sind hier dadurch möglich, daß aus den Signalen 
der einzelnen Strahlen des Gamma-Densitometers zunächst die 
Strömungsform bestimmt wird und anschließend die Querschnitts- 
mittelung des Dampfvolumeranteils in Abhängigkeit von der 
Strömungsform vorgenommen wird /24/. 


Abb. 9 zeigt die entsprechenden Ergebnisse gemessen mit dem 
LOFT-MeBeinsatz: Die Streuung der Meßwerte ist erheblich. Ab- 
weichungen bis zu einem Faktor zwei und mehr treten afif. Es 
ergeben sich für die einzelnen Meßkombinationen folgende 
Mittelwerte x und Standardabweichungen c : 


x G 
br Ber 1,36 0.43 


G -p0 Ref 0,95 0,20 


Gro G r 073 0,32 


Von den drei Kombinationsmöglichkeiten ergibt die Kombination 
Gamma-Densitometer-Stauscheibe die besten Ergıhnisse. 


In Abb. 9 sind keine MeBpunkte vorhanden, bei denen ein MeB- 
wert "falsch" gemessen wurde aufgrund von z.8. Verstárker- 
drift oder Messen auderhalb des Gerätemeßbereichs. Die Ab- 
weichungen sind dadurch bedingt, daß zum einen die ver sende- 
ten Beziehungen nur für homogene Strómungsform gültig sind, 
zum anderen die Stauscheibe und der Turbinenradzähler in 
geuissem Maße örtliche Werte anzeigt. 
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Die Abweichungen, die durch die Meßwertbez 
hervorgerufen werden, sind in Abb. 10 darg 
Hierfür wurden in den Gleichungen (7)-(9) 

ሃን die Gleichung (4)-(6) eingesetzt und di 
allgemein gültige Beziehung für die Massen 
dividiert. Abb. 10 zeigt die rechnerisch b 
stromdichte-Verháltnisse als Funktion des ` 
háltnis aus Phasengeschwindigkeit des Gase: 
digkeit der Flüssigkeit). Parameter ist de 


a 


Die Auswertebeziehungen, die die MeBwerte ( 
meters und Turbinenradzählers sowie der St; 
Turbinenradzählers kombinieren zeigen eine 
vom Schlupf. Dagegen ist die Auswertebeziet 
nation Gamma-Densitometer-Stauscheibe verg | 
Schlupf abhängig. Diese Tatsache erklärt zu 
Ergebnis für diese «ombination. 


Wie erwáhnt, überstreichen die hier betrach 
Gebiet derjenigen Strómungsformen, die dure 
ausgeprágte Dichteschichtungen besitzen. Da 
Strömungsgröðen in einem Rohrquerschnitt im 
der Vertikalen abhánoen. 


Mit den Gamma-Densitometor wird ein verglei 
wert der Dichte bzw. des Dampfvolumenantei] 
Menrstrahl-Gamma-Densitometer kann darüberh 
über die Dichteverteilung im horizontalen R 
die hier diskutierten Versuche war es sehr - 
Hilfe des Gamma-Densitometers die Hóhe der | 
bestimmen, die bei einer idealen Schichtens: 
ware (Abb. 11 unten). Es zeigte Sich, daB d 
trennfláche (ሃ/4)፲፥ geeignet ist, um einen: 
hang zwischen lokalen und integralen Ne iwer: 
und des Turbinenradzählers herzustellen / 15; 





iehungen (7) bis (9) 
estellt (aus /15/). 
für 0, ( ov? lan und 

e Werte durch die 
stromdichte G1.(10) 
estimmten Massen- 
Schlupfes S (Ver- 

5 zu Phasengeschwin- 
P Dampfvolumenantei | 


les Gamma-Densito- 
luscheibe und des 
starke Abhängigkeit 
jung für die Kombi- 
eichsweise wenia vom 
im Teil das giinstigere 


teten Versuche das 

h Schwerkrafteinflu3 
s bedeutet, dap die 
wesentlichen von 


chsweise guter Mittel- 
5 bestimmt. čin 

inaus noch Information 
phr liefern /24/. Fiir 
vorteilhaft, mit 
dasseroberfläche zu 
trömung vorhanden 

lese Höhe der Phasen- 
sinfachen Zusammen- 


ten der Stauscheibe 
f 








Wasserdampf (Vol. -“/o) 


Dieser Zusammenhang wird auch verwendet, um die LOFT-Massenstrom- 
instrumentierung zu kalibrieren. In Abb. 11 ist dies für die 
Kombination Stauscheibe-Gamma-Densitometer dargestellt: Die Be- 
ziehung (8) wird mit einem Koeffizienten C multipliziert, der in 
einfacher Weise von der Höhe der Phasenstromfläche abhängt. 

Die auf diese Weise bestimmte Massenstromdichte ist in Abb.12 
dargestellt. Die Streuung der Meßpunkte wurde deutlich reduziert; 
die mittlere Abweichung des Massenstromdichteverhältnisses 
beträgt # 4 X, die Standardabweichung ፦15 Z. 











In der Abb. 13 sind die Ergebnisse der PHEBUS-Instrumentierung 
aufgetragen. Zusätzlich zu der Variation von Massenstrom, Dampf- 
gehalt und Druck wurde hier auch die Strömungsrichtung variiert. 
Die Genauigkeit der Resultate für nach oben gerichtete Strömung 
ist gut, die Abweichung des Mittelwerts beträgt ca. 1 %, die 
Standardabweichung 6 %. Für abwärtsgerichtete Strömung treten 
größere Meßfehler auf. Hier ist die mittlere Abweichung - 16 %, 
die Stand ırdabweichung 12 %. 


Die Abweichungen werden dadurch verursacht, daß die Annahme 
gleicher Dichten in den Meßquerschnitten nicht genau zutrifft. 
Aufgrund der Beeinf!ussung der Phasenverteilung durch Strömungs- 
umlenkungen im Rohrsystem und durch die Meßgeräte selbst, kann 
der Dampfvolumenanteil und damit die Dichte in den Meßauer- 
Schnitten verschieden sein. Dieser Sachverhalt läßt sich auch 
durch den Schlupf S erfassen, der sowohl für die Venturidüse, 
als auch für den Turbinenzähler mit Hilfe der Peferenzwerte 

der Kreislaufinstrumentierung berechnet wurde /25/. Es zeigte 
sich, daß der Schlupf für jedes Meßinstrument nur schwach vom 
Massenstrom und Dampfgehalt abhängig war, jedoch eine deutliche 
Abhängigkeit vom Druck und der Strömungsrichtung besaß. Abb. 14 
zeigt die für jeden untersuchten Druck gemittelten Werte, 
aufgetragen über dem Dichteverhältnis der Phasen bzw. dem Druck. 













Das Prinzip der Kalibrierung dieser Instrumentierung beruht 
nun darauf, diese Abhängigkeit des Schlunfes in den entsprechen- 
den Auswertebeziehungen zur Bestimmung der Dichter zu verwenden. 





die Werte für abwärtsgerichtete Strömung, die erhebliche Ab- 
weichungen besaßen, wurden deutlich verbessert (Abweichung des 
Mittelwerts 3 %, Standardabweichung 13 %). 


6. Besondere Merkmale der Meßgeräte 





Im folgenden werden stichwortartig die Ergebnisse zusammengefaßt 
sowie sehr vereinfachend andere typische Merkmale der Meßgeräte 
kurz angesprochen: 


TMFM: 

Die Vorteile bestehen in der einfachen Bestimmung des Massenstroms 
da nur eine Meßgröße benötigt wird und der Unabhängigkeit von der 
Phasenverteilung (Strómungsform). 

Als dachteil ist zu werten, daB die Strómung stark gestórt wird 
und das MeBgerát einen Druckverlust verursacht, der bis zu einigen 
bar betragen kann. Der mechanische Aufbau des TMFM ist sehr auf- 
wendig. 





Radionukiid-Verfahren 





Als Vorteil ist zu sehen, daß die Phasengeschwindigkeiten direkt 
gemessen werden und damit auch der Dampfanteil der Strömung be- 
stimmt werden kann, sowie die weitgehende Unabhängigkeit des Meß- 
verfahrens von der Phasenverteilung. 

Ein Nachteil besteht darin, daß die Verwendung von radioaktiven 
Isotopen einen erheblichen Aufwand mit sich oringt (Isotopen- 
Produktion - Transport, Strahlenschutz). Die Anforderungen an die 
MeBtechnik und Signalverarbeitung sind hoch. 


Stauscheibe und Turbinenradzähler (DTT) und Gamma-Densi tometer 





Die Abb. 15 zeigt die Ergebnisse mit diescr Kalibrierungsmethode: 





Fin Vorteil ist, daB der DTT auch ir erschwert zugánglichen 
Geometrien eingesetzt werden kann (Downcomer, Gitterplatte) 


Nachteilig ist, daB eine starke Abhängiakeit der Meßwerte vom 
Strömungszustand vorhanden ist, was eine aufwendige Kalibrierung 
notwendig macht. 


Venturidüse und Turbinenradzähler 





Der Vorteil liegt hier im einfachen Aufbau und der einfachen Me&- 
wertverarbeitung. 
Nachteilig ist auch hier die Notwendigkeit der Kalibrierung. 


7. Kalibrierung der MeBverfahren in instationärer Strömung 





Die bisherigen Ausführungen beschränken sich auf Untersuchungen 
in statiönärer Strömung obwohl die Meßverfahren für Experimente 
entwickelt wurden, in denen sich der Strömungszustand zeitlich 
stark ändern kann. Die Begründung für diese Vorgehensweise war, 
dai 
a) eine Referenzmessung des Massenstroms (z.B. durch Messung 
von einphasigen Massenstrómen stromauf de” Mischkammer, 
wie in Kap. 4 beschrieben) nicht eristierte. 


b) eine Analyse der Abhängigkeit der MeBsignale von bestimmten 
Phasenverteilungen (Strómungsformen) nur múglich ist, wenn 
diese Phasenverteilungen über einen längeren Zeitraum (stationär) 
eingestellt werden und deshalb gemessen werden können (z.B. 
mit Impedanzsonde und Gamma-Densitometer). 


Auf diese Weise wurden Erkenntnisse gewonnen, die zum Teil 

direkt auf instationäre Strömung übertragen werden können. So ist 
zu erwarten, daß ein Meßverfahren, dessen Meßgenauigkeit sich in 
stationärer Tests a's prakt"sch unabhängig von der Phasenvertei- 
lung erwiesen hat, auch bei instationären Tests keine merkliche 
Abhängigkeit von anderen Pnasenverteilungen zeigen wird. Auf der 
anderen Seite ist diese Obertragbarkeit nicht a priori zulässig 
bei Meßverfahren, deren Genauigkeit vom Strömungszustand abhängt. 
Auch für die erstgenannten Meßverfahren gibt es Gründe, die Tests 
in instationärer Strömung notwendia machen. Hier ist zu erwähnen 
die starke mechanische Beanspruchung der Meßgeräte bei Beginn 
eines Blowdowns, der Einfluß starker Temperaturtransienten auf die 
MeBgenauigkeit und die zeitliche Auflösung der Meßverfahren. 


Aus diesen Gründen wurde der Aufbau eines Versuchsstandes für 
instationäre Strömung (Blowdownversuchsstand) in Angriff ge- 
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nommen /26/, der 1981 für Tests verfügbar sein wird. Die Referenz- 
messung soll dabei mit einem TMFM erfolgen, das speziell für 

diese Betriebsbedingungen ausgelegt wurde. Als zusätzliche 
Instrumentierung in der Teststrecke werden Gamma-Densitometer 
eingesetzt, die sich in stationärer Strömung bewährt haben 

und für instationäre Strömung modifiziert werden. 
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Tests of Two-Phase Mass Flow Rate Instrumentation 


Tab. 1 
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True Mass Flow Meter (TMFM) (KIK-IRE-IT) 














Prinzip einer radialen Strómungsmaschine (schematisch! 
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FESTLEGUNG DER EICHFAKTOREN 
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(y/d)= =0.5 (1-cos lei 2) ) 














Abb.11: Abhängigkeit des Kalidrierungskoeffizienten 
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Zusammenfassung 











Es wird eine Obersicht gegeben über einige im Institut für 
Reaktorentwicklung bereitgestellte Rechenverfahren zur 
Fluid- und Strukturdynamik. Folgende Anwendungsfälle werden 











beschrieben: 


- Belastungen der Reaktordruckbehälter-Einbauten (insbeson- 
dere Kernmantel) eines DWR beim Bruch einer Kühlmittel- 























leitung .m kalten Strang. 


- Belastungen der Core-Stützen(Steuerstabführungs rohre) 
eines DWR beim Bruch im heißen Strang. 






Belastungen von kugelfórmigen Containments mit Druckunter- 
drückungssystem (Paulinie 69) durch Dampfkondensation. 








Schwingungen von Zylinder- und Kugelschalen. 





Kompressible und Zweiphasen-Strómung in komplizierten 






Geometrien. 
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Auf die Bedeutung einer verbesserten Analyse mit derartigen, 
relativ aufwendingen Rechenverfahren wird hingewiesen. 
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Unsere bisherigen Arbeiten / ;, 2, 3 / konze 
ausschließlich auf Blowdown-Störfälle. Es wu 
ang>nommen, daß die Hauptkühlmittel-Leitung | 
zwar 


- daß der Bruch an der ungünstigsten Stelle 
nahe am Einlaßstutzen oder Auslaßstutzen : 
behälters), 


- daß der ungünstigste Bruchzeit,erlauf vor 
artiger Bruch), 


- daß die durch den Bruch entstehende Ausst 
vollen Rohrquerschnitt entspricht. 


Man sollte sich dies vor Ausen halten. Ist di 
Hauptkühlmittel-Leitung schon sehr unwahrsch: 
ja alles getan, um ihn aus*uschlieBen -, so | 
der ungünstigsten Stelle, darüber hinaus noc! 
Bruch mit voller Querschnittsfreisetzung, no: 
ordnungen unwahrscheinlicher. 


Trotzdem, für diese extremen Ereignisse, gen: 
dieser Ereignisse, wurden vom Institut für Re 
Methoden entwickelt, die zumindest in einzeli 
die derzeit üblichen Verfahren hinausgehen. | 
einmal um Methoden zur Fluiddynamik. Es inter: 
sondere die Druckverteilung an den Fluid-Ran: 
stimmt die Strukturbelastung. Strömungszustä: 
sind dagegen von geringem Interesse. Es hand 
Methoden zur Strukturdynamik. Hier interessit 
den Verformungen. Sind diese Verrcrmungen to! 
sie zu Brüchen? ist insbesondere die Nachwärı 
Abschaltbarkeit des Reaktors gewührleistet? 








ntrierten sich fast 
rde bekannterweise 
brechen möge; und 


erfolgt (meist 
des Reaktordruck- 


liegt (meist schlag- 


römfläche dem 


sr Bruch an einer 
:inlich - man hat 
ist ein Bruch an 
ı ein schlagartiger 
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nen Punkten über 

zs handelt sich 
:ssiert hier insbe- 
iflächen. Sie be- 
ide im Feldinneren 
211 sich ferner um 
ren die auftreten- 
lerierbar? Führen 
peabfuhr und die 





In vielen Füllen kónnen das fluiddynamische und struxturdynami - 
sche Teilproblem nicht unabhängig voneinander gelöst werden. Bei 
der Behandlung des fluiddynamischen Teilproblems müssen cie 
Wandverschiebungen bzw. -bewegungen bekannt sein. Sie kängen aber 
vom strukturdynamischen Verhalter ab. Dieses wiederum kann nur 
beschrieben werden, wenn die Beiastungen aus der Fluiddynamik 
bekannt sind. "Die Katze beißt sich also in den Schwanz." Mit 
anderen Worten, es handelt sich hier um ein Problem der Fluid- 
Struktur-Wechselwirkung. Es kann nur als ein einziges, zusammen- 


hängendes Problem gelöst werden. 


1. Computerprogramm FLUX [47 





Diese Methode erlaubt die Berechnung der Belastungen und Bean- 
spruchungen von Reaktordruckbehälter-Einbauten beim Blowdown 
durch den Einlaßstutzen. Die Anwendung erfolgt für den Druck- 
wasserreaktor. Zuerst eine kurze Charakterisierung des Modelles: 
FLUX liefert eine dreidimensisnale Beschreibung des Problems. 
Die Fluid-Struktur-Wechselwirkung mit dem elastischen Kernmantel 
wird berücksichtigt (das Schalenmodeil für den Kernmante! wird 
später erläutert). Die Fluid-Struktur-Wechselwirkung mit einem 
homogenisierten Coremodell und einem nachgiebigen [.eaktordruck- 
behälter ist in Vorbereitung. (Zur Zeit wird das Core lediglich 
durch eine Punktmasse am unteren Ende des Kernmantels simuliert 
und der Reaktordruckbehälter wird als starrer, eingespannter 
Körper behsudelt.) Die Blowdown-Strömung wird als kompressible 
Potentialströmung behandelt. Unter diesen Bedingungen - drei- 
dimensionales Modell, Fluid-Struktur-Wechselwirkung, kompressib- 
le Strömung - führt eine ausreichend feine Diskretisierung zu 
einem sehr hohen Rechenaufwand. Dieser ist nur dann tragbar, 
wenn ausgefeilte numerische Techniken eingesetzt werden. Die 
Entwicklung soicher Techniken war das Hauptproblem bei der Er- 
stellung des Computerprogrammes FLUX. Was sind nun die Merkmale 
dieser Techniken? Die Bewegungsgleichungen für Fluid- und 
Strukturdynamik werden simultan und imp2izit integriert. Simultan 
bedeutet keine Iteration an der Fluid-Struktur-Kontaktflache. 








Implizit bedeutet keinen Zwang zu sehr kleinen Zeitschritter.. 

Es bleibt, pro Zeitschritt ein freilich sehr großes, lineares 
Gleichungssystem zu lösen. Dazu wird ein sogenannter schneller 
elliptischer Lóser eingesetzt. Er zieht Vorteile aus der be- 
sonderen Struktur der Koeffizienten-Matrizen der linearen 
Gleichungssysteme. Diese besondere Struktur hàng- von der Geo- 
metrie des Problems ab. Dies bedeutet, daß beim Finsatz schneller 
elliptischer Löser in besonderer Weise auf die Geometrie des 
Problems Rücksicht genommen werden muß. 


Nun einige Bilder zur Demonstration: Bild ! zeigt ein Isobaren- 
feld, wie es sich kurz nach Biowdown-Auslósung einstellt. Man 



















Bild 1 


Berechnetes Isobarenfeld zur Zeit 

t = 20 ms in einem Druckwasserreaktor. 
Der Abstand zwischen den Isobaren be- 
trágt 1 bar. Der Anfangsdruck ist 

Po * !57 bar. Die gestrichelten Iso- 





Aas dem Verlauf der Isobaren im Innen- 
bereich erkennt man einen Einflu& der 
Ke rnmante lbewegung. 











erkennt deutlich die nichtrotationssymmetrische Druckverteilung. 
Bemerkenswert ist auch die Geometrie: diinner Kernmantel, unter- 
schiedlich groBe Dimensionen des Fluidbereiches (kleiner Stutzen- 
durchmesser, kleine Ringraumbreite, groBe Hóhe). Eine wesentliche 
Voraussetzung für die Zuverlüssigkeit solcher Rechnungen ist eine 
ausreichend feine Diskretisierung. Dies geht aus Bild 2 hervor. 
Es zeigt die maximale Kernmantelspannung als Funktion der Zeit. 






baren gehóren zu Drücken unter 150 bar. 














Bild 2 Maximale Vergleichsspannung als Funktion der Zeit für 
verschiedene Diskretisierungsfeinheit N1 bis N4. Die Zah- 
len der Fluidmaschen/Struktureigenschwingungen sind 


N1: 170/15 N2: 1022/88 N3: 4326/271 N4: 12663/596. 























Das Maximum tritt nic:. immer an derselben Stelle des Kernmantels 
auf. Die vier verschiedenen Kurven wurden mit unterschiedlichem 
Diskretisierungsaufwand erzielt: Kurve N1 geringer Aufwand, N2 
hóherer Aufwand usw. Wir erkennen, daB erst der Diskretisierungs- 
aufwand bei Kurve N3 ausreichend ist. Es wurden hier etwa 4000 
Fluidmaschen und 270 Eigenschwingungs formen des Kernmantels ver- 
wendet. Die Rechenzeit auf der IBM 370/3033 betrügt heute etwa 
eine halbe Stunde. Bild 3 zeigt einen Vergleich zwischen einer 
Drw.kmessung bei Blowdown-Versuchen am HDR und Rechnungen mit 
FLUX. Einmal wurde in FLUX die Fluid-Struktur-Wechselwirkung be- 
rücksichtigt und ein anderes Mal nicht. Man erkennt sofort, ohne 
Berücksichtigung der Fluid-Struktur-Wechseiwirkung ist der berech- 











nete Zeitverlauf vóllig falsch. 




















Bild 3 
Druckdifferenz am Kerr- 
mantel im HDR-Blowdown- 
Experiment V31 













Druckdifferenz Í Þar) 





























2. Computerprogramme SING! und SING-S /^5 7 











Diese Codes werden eingesetzt für trarsiente Fluiddynamik- 
Probleme. Sie sind nicht zugeschniiten auf spezielle Probleme 
wie beispielsweise FLUX. Weitgehend beliebige, dreidimensionale 
Geometrien kónnen behandelt werden. Fluid-Struktur-Wechselwir- 
kung kann berücksichtigt werden. In diesem Falle müssen die 
Steifigkeitsmatrizen der Strukturen vorgegeben werden. Sie kónnen 
beispielsweise mit Hilfe der bekannten Finite-Elemente-Me thoden 
bestimmt werden. Selbstverständlich hat die Anwendbarkeit auf 
weitgehend beliebige, dreidimensionale Geometrien auch ihren 
Preis: Die Strömung muß als inkompressitle Potentialstrómung 
behandelt werden, d.h. Welienausbreitungen erfolgen nicht mit 

der endlichen Schallgeschwindigkeit, sondern unendlich schnell. 
Dies ist akzeptabel, wenn die Wellenausbreitungszeiten ohnehin 
klein sind gegenüber den Problemzeiten. Bei inkompressibler 
Potentialstrómung und Fluid-Struktur-Wechselwirkung treten übri- 
gens nur die Struktur-Schwingungsfreiheitsgrade auf. Insbeson- 
dere bei Annahme eines inkompressiblen Fluids lassen sich die 
sogenannten Randintegral-Verfahren - auch bekannt unter dem 

Namen Singularitäten-Verfahren oder Panel-Verfahren - sehr 
effektiv zur Lósung einsetzen. Das dreidimensionale Fluiddynamik- 
Problem geht hierbei in ein zweidimensionales Probiem in den 
Fluid-Randflächen über. Es ist deshalb nur eine zweidimensionale 
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Diskretisierung der Fluid-Randflächen notwendig. Bei den uns 
interessierenden Ergebnissen kommt es aber gerade auf die Fluid- 
Randflächen an. Nur die dort auftretenden Drücke und Normalver- 
schiebungen sind von Interesse. Die Codes SING! und SING-S 

sind relativ leicht anzuwenden und der Lésungsalgorithmus ist 
relativ gut zu durchschauen. Die Codes wurden bisher für zwei 
Probleme aus der Reaktorsicherheit eingesetzt. 


Randflächen des Fiuidbereiches wurden mit den hier — ^ 
Panels diskretisiert. Eine relativ große Zahl der Corestützen 
wurde berücksichtigt. Sie sind in Bild 4 engedeutet. Selbstver- 
stindlich stellen auch die Stützen Fluid-Randfláchen rd Ihre 
Diskretisierung mit Hilfe von Panels wird in Bild 4 ar. ge- 
zeigt. Nach ersten Rechnungen liegen die Spannungen w den am 
höchsten beanspruchten Corestützen im Bereich der €— 
Einige der gemachten Annahmen nüssen allerdings noch überprüft 


— werden. Die Rechenzeit betrug etwa 30 Minuten. 
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Ziel ist die Beanspruchung der Stützen. Bild 4 zeigt eine verein- ; 1 


fachte Geometrie des oberen Plenums mit den Stützen. Die äußeren 


Die Anregung erfolgt durch Dampfkondensation ım Wasserpool des 
Druckunterdrückungssystems. Ziel ist die w "we 
fórmigen Containmentschale. Bild 5 zeigt su. 8 "ሸዊ | 
des Containments. Die Fluid-Randfláche wurde mit as me vu 
diskretisiert (das Verhalten der Kugelschale wurde w መ 
spezieller Modelle beschrieben. Diese werden w. näher er- 
láutert). Die berechneten Eigenfrequenzen der Kugelschale, "- 
Berücksichtigung der Fluidtrágheit des Wasserpools, beginnen mit 
SO Hz und liegen dann dicht an dicht. Die Figenfrequenzen m:t 


Serücksichtigung der Trágheit des Wasserpcols beginnen mit Werten 





ini | s en auch bei Messungen in 
um 10 Hz. Einige dieser Frequenzen wurd 









der Anlage Brunsbüttel festgestellt. Der Einflu& der Fluid- 






Struktur-Wechselwirkung auf das Schwingungsverhalten ist also 






sehr groß. Als Ergebnis dieser Untersuchungen wurden nicht wie 






üblich Containment-Spannungen ermittelt. Dazu wáre námlich die 






| ignisse notwendi ewesen. 
Vorgabe maximaler Kondensationsereignis g g 






Dies war aber wegen des statistischen Charakters dieser Ereig- 





nisse nicht móglich. Es wurden vielmehr die eben noch zuldssigen 
> = 






| | f _ Leg bag 
Kondensationsereignisse - auch Kombinationen davon ermittel 






Ü ; en gerade aus. 
| | Sie schópfen die zugelassenen Spannung g 
Bild 4 Oberes Plenum eines Druckwasserreaktors. 








Die äußeren Wände sind angenáhert durch recheckige 
Konturen. 


Wánde und Stützen sind diskretisi 


eines Randintegral-Verfahrens. 


ert zur Árwendung 












Bild 2 


Kugelfórmiger Sicherheitsbehälter 
— — Wetten L. |] I} mit Druckunterdrückungssystem 
Tube unter Blrædown-Belastung. 
Die Oberfláche des Wasserpools 
A des Druckunterdrückungssystems 

i | ist diskretisiert zur Anwen- 
dung eines Randintegral-Verfah- 
rens. 
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Dezu wird das in Bild 6 dargestellte einfache Beispiel verwendet. 











Es handelt sich um ein dickwandiges Rohr, gefüllt mit einem 

Fluid. An der einen Stirnseite ist der Druck p aufgeprägt, an 
der anderen Seite soll plótzlich eine Druckabsenkung erfolgen. 
In der Mitte des Fluids befindet sich ein starrer Stab (Fluid- 


Struktur-Wechselwirkungsphänomene spielen also keine Rolle). 


Was passiert nun bei der Druckabsenkung? Zwischen den beiden 










Bild 6 
Zur Belastung eines 
starren Stabes in 
einem transienten, 
Mittelwert Schnligeschmmhgie kompressiblen und | 
Last - fesk | K  inkompressiblen Fluid- 
x f feld 
















































Stirnflächen laufen Druckwellen mit der Schallgeschwindigkeitt 
Passiert eine solche Druckwelie den Bereich des 











hin und her. 
Stabes, so ruft sie einen Lastpeak hervor, der zwar sehr kurz, 


dafür aber recht hoch ist. Tatsächlich kommt es 2lieine auf die 












Hóhe dieser Lastpeaks gar nicht an, wenn nur di: Wellenausbrei- 






tungszeiten deutlich kleiner sind als andere, für das Problem 








charakteristische Schwingungszeiten. Unter dieser Voraussetzung 





ist allein der Belastungsmittelwert ausschlaggebend. 





Denkt man sich nun die Schallgeschwindigkeit erhöht, beispiels- 





weise auf das Dreifache, denn sind die Lastneaks nur ein Drittel 







so breit und die Abstünd^ zwischen diesen Lastpeaks auch ein 
Drittel so lang wie zuvor. Der Lastmittelwert aber bleibt derselbe. 
Des heißt, sind die Wellenausbreitungszeiten nur deutlich kleiner 


als die sonstigen charakteristischen Schwingungsze:ten, kommt es 














Bild 7 zeigt einige berechnete Eigenschwingungsformen für die 


Zylinderschale. Man beachte die starke, auf engen Raum konzentrier- 


also nur auf den Lastmittelwert an, so kann die Schall geschwindig- 
keit beliebig erhëht werden. Insbesondere kann sie auch unendlich 
groß sein. Das bedeutet aber, dic Fluidkompressibilitát kann ver- 
nachlässigt werden. 







Im Falle der Blowdewn-Belastung des Kernmantels wurde dies anhand 
von Rechnungen mit FLUX bestätigt. Ohne Berücksichtigung der Kom- 
pressibilität ergaben sich nur wenig höhere Beanspruchungen als 

mit Berücksichtigung der Kompressibilität / 4 7. 







Die bisher besprochenen Programme FLUX und SING sind für größere 
Probleme gedacht. Daneben wurden aber auch fortgeschrittene Metho- 
den für Teilprobleme entwickelt. 


3. Computerprogramme CYLDY3 / 7 7, SPHERE-DYNA / 8 7 





Diese Programme werden für schalendvnamische Probleme eingesetzt. 
CYLDY3 beschreibt das dyramische Verhalten dünner Kreiszvlinder- 
schalen. Beide Schalenenden sind ringversteift. Das eine Schalen- 
ende kann entweder starr oder elastisch gebettet scin, das ande- 
re ist fr:i beweglich. CYLDY3 wurde bisher zur Untersuchung der 
Kernmantelschwingungen beim Blowdown eingesetzt. SPHERE-DYNA 
beschreibt das dynamische Verhalten dünner Kugelschalen. An zwei 
planparailelen Kreisen sind spezielle Randbedingungen vorzugeben. 
SPHERE-DYNA wurde bisher fiir Untersuchungen von Containment- 
Schwingungen, angeregt durch Kondensationsvorgánge ir. Druckunter- 
drückungssystem, eingesetzt In beiden.Codes, CYLDY3 und SPHERE- 
DYNA, ist es gelungen. fiir die Eigenschwingungen analytische 
Lósungen herzuleiten. Fs wurden hierbei die ungekúrzten Schalen- 
differentialgleichungen von Flügge zugrunde gelegt. Da analytische 
Lósungen nicht an Maschensysteme gebunden sind, kónnen mit CYLDY3 
und SPHERE-DYNA sehr hohe órtliche Auflósungen erzielt werden. 















te Biegung an der Einspannstelie. Bei der Anwendung klassischer, 
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Bild 7 

Radiale ( ) und 
azimutale (----) Ver- 
schiebungen für die 








erste Biege-Eigenschwin- 
gungsform 











ösungsverfahren würde eine solche Beschreibung einen 
itunter das Mögliche übersteigenden Diskretisierungs- 
dern. Biid 8 zeigt die für die Zylinderschale berech- 
equenzen. Sie sind aufgetragen über der Umfangsord- 
ingungen. Eigenfrequenzen, die durch Linien mitein- 
en sind, gehören zu ähnlichen axialen Schwingungs- 
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Bild 8 

Eigenfrequenzen des HDR- 
Kernmantels, aufgetragen 
über der Umfangsordnung. 
Die Linien kennzeichnen 
Lësungen gleicher Axial- 
ordnung. 
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formen. Entsprechende Kurven wurden auch fúr die Kugelschale be- 

rechnet. Dort liegen die Frequenzen noch wesentlich dichter bei- 

einander. Bei den gewáhlten Abmessungen betrágt die kleinste Fre- 
quenz etwa 50 Hz. Im Bereich bis 70 Hz liegen bereits etwa 

300 Eigenfrequenzen. Bei dynamischen Analysen benutzten wir bis 


zu 2000 Eigenfrequenzen. 






Bei dünnen Kugelschalen stellt msn übrigens ein überraschendes 
"lokales" Verhalten fest: Eine in einem kleinen Bereich angeregte 
Kugelschale macht wesentliche Schwingungsausschláge nur in dem 
selben Bereich. Die Nachbarbereiche bleiben, wie in Bild 9 er- 
sichtlich, relativ ruhig. Diesen Effekt erhált man nur, wenn sehr 





Bild 9 
Dynamische Reaktion 
einer lokal belasteten 


Kugelschale 











viele Eigenfrequenzen und Eigenschwingungsformen ` rücksichtigt 
werden. Bei der Verwendung von Standard-Finite-Elemente-Me thoden 
ist dies wegen des beschrünkten Diskretisierungsaufwandes nicht 


móglich. 





CYLDYZ wurde getestet durch Vergleich mit anderen Codes. In einem 
Falle war die Abweichung in den Eigenfrequenzen weit unterhalb 
von einem Prozent. SPHERE-DYNA wurde getestet durch Vergleich mit 


Messungen am Reaktor Brunsbüttel. " 





4. Computerprogramme FLUST / 9 7, DRIX-2D / 10_7 











Diese Codes werden eingesetzt für Teilprobleme aus der Fluid- 
Dynamik. FLUST gestattet die Beschreibung zweidimensionaler 
kompressibler Strómungen in beliebig areinander gekoppelten, 
rechteckigen Bereichen. Dies ist ein Fortschritt gegenüber den 
meisten bisher üblichen Verfahren, bei denen die Strömungsfelder 
durch Rohrleitungsnetze simuliert wurden. In diesem Falle handelt 
es sich um eindimensionale oder null-dimensionale, aneinander ge- 
koppelte Bereiche. Bild 1O zeigt solche zweidimensionale, anein- 
ander gekoppelte Bereiche für Reaktor-Geometrie. K ist der abge- 
wickelte Ringraum (das ist der Fluid-Bereich zwischen Reaktor- 


































































QUERSCHNITT OURCH FLUST - 
DEN HDR - RDB MODEL! GEOMETRIE 

















Bild 10 Beschreibung der Fluid-Bereiche im HDR-Druckbehälter 
durch aneinander gekoppelte, zweidimensionale Bereiche 









Druckbehálter und Kernmantel), S ist der Stutzenbereich und C ist 
der äußere Bereich im Containment. Beim Blowdown strömt das Fluid 
vom Ringraum R durch den Stutzen S in das Containment C. DRIX-2D 
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beschreibt transiente Zweiphasen-Strómungen im thermodynamischen 
Nicht-Gleichgewicht mit Schlupf zwischen den Phasen Die Ge- 
schwindigkeitsfelder der Flüssigkeit und des Dampfes werden g^- 
trennt voneinander behandelt. Dies geht aus Bil; 11 hervor. Darge- 
stellt ist die Strómung in einem T-Stück. Wasser und Dampf haben 
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Bild 11 Strëmungsfeld einer Zweiphasenstrómung mit Schlupf in 


einem T- túck 


unterschiedliche Strómungsgeschwindigkeiten, vor allem aber unter- 
schiedliche Strömungsrichtungen. Es liegt also ein Schiupf zwischen 
beiden Phasen vor. Der Dampf läßt sich leichter umlenken als das 
Wasser. Zur Zeit wird eine Kopplung des Zweiphasen-Codes DRIX-2D 
mit FLUX vorbereitet. DRIX-2D wird dabei zur Berechnung der Stró- 
mung mit Blowdown-Stutzen, FLUX für die Berechnung der Strómung 

und Fluid-Struktur-Wechselwirkung im Druckbehilter eingesetzt. 


Generelle Bemerkungen 





Computercodes können sehr leistungsfähige Werkzeuge sein, wenn 
sie sachgemäß angewandt werden. Insbesondere sollte man die ein: 
schränkenden Annahmen niemals aus dem Blickfela verlieren. 
Effekte, die man im Prinzip schon vorher kannte, werden durch die 
Codes quantitativ erfa8t. Bei Effekten, die man dagegen vorher 











nicht kannte, ist äußerste Vorsicht geboten. Gt 
gungen sind unbedingt erforderlich. Die gefundi 
sollte man erst dann als echt akzeptieren, weni 
liegenden Phänomene qualitativ durchschaut sin 


Auf welche Art sollen verbesserte Rechenmethodt 
bei der Stórfallanalyse verbessern? Was nutzt « 
wenn man ein Teilproblem mit verbesserten Meth: 
zent genau beschreiben kann, auf der anderen Se 
gangsparameter wesentlich größere Unsicherheit: 


Man halte sich dazu die bisherige Situation voi 
oberen Teil des Bildes 12 charakterisiert. Im : 
heitsanalyse von der Definicion der Eingangspaı 
Schlußaussagen - dies ist die Abszisse in Bild 
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Unsicherheiten, Unschärfen und Fehler in allen 
stadien auf. Die Fehler sind über die ganze Sic 
mehr oder weniger stark verschmiert. Zur Gesamt 
Genauigkeit benótigt man Wissen auf allen mógli 
chen gleichzeitig. Man muß die Sicherheitsphilc 








:nauere Oberle- 


nen Effekte 
1 die zugrunde 


l. 


^n den Wissensstand 
'S beispielsweise, 
den auf ein Pro- 
ite aber die Ein- 
"n aufweisen? 


r Augen. Sie ist im 
luge eiuer Sicher- 
rameter bis zu den 
12 - treten 


| 12 

Erlüuterung der Ent- 
lung von Ungenaui g- 
ten währen, des Ab- 
es einer Sicher- 
‚sanalyse 


möglichen Zwischen- 
the rheitsanalyse 
beurteilung der 
chen Teilberei- 
sophie kennen, hin- 


zu kommen praktische Erfahrungen, aber auch Kenntnisse der mathe- 
matischen und numerischen Methoden, schließlich auch Kenntnisse 
der Werkstoffe. Dies ist für eine zuverlässige Beurteilung sehr 
hinderlich. 


Die Situation ist anders, wenn ausreichend genaue Methode für 
Teilbereiche zur Verfügung stehen. Die im Prinzip unvermeidli- 
chen Fehler und Unschürfen werden an strategisch günstigen Stel- 
len konzentriert. Dazu ein Beispiel: Containment mit Druckunter- 
drückungssystem, Anregung durch Kondensatio*svorgünge im Wasser- 
pool. Die gekoppelten fluid-strukturdynamischen Vorgänge konnen 
mit Hilfe der von uns entwickelten Codes beschrieben werden / 6 /. 
Es wurden nun nicht wie sonst üblich die berechneten mit den zu- 
lássigen Spannungen verglichen, denn bei einem solchen Vergleich 
liegt das eigentliche Problem nicht. Es wurden vielmehr unter Zu- 
grundelegung der zulássigen Spannungen zulássige Kondensations- 
ereignisse bestimmt. Diese sind nun mit den tatsáchlich beobach- 
teten Kondensationsereignissen zu vergleichen. Dies ist sinnvoll, 
weil die Intensitát der Kondensationsereignisse jene physikali- 
sche Größe ist, die naturbedingt mit den größten Unschürfen behaf- 
tet ist. Wer das Gesamtprobiem beurteilen will, braucht sich nun 
nicht mehr mit den gekoppelten fluid-strukturdynamischen Problemen 
zu befassen, sondern er kann sich auf das eigentliche Problem der 
Dampfkondensation konzentrieren. 


Außerdem, Fehler ist nicht gleich Fehler. Ungenauigkeiten und 
Vernachlüssigungen bei der Lósung von Differentialgleichungen sind 
schwer zu überschauen. Variationen bei den Randbedingungen oder 
sonstigen Eingabeparametern sind dagegen physikalisch interpretier- 
bar. Die Konsequenzen solcher Unschürfen kann der Ingenieur im 
allgemeinen verstehen. Deshalb ist es durchaus sinnvoll, bei mathe- 
matischen Prozeduren hohe Genauigkeiten zu verlangen. Dies haben 
wir insbesondere durch die Entwicklung der schalendynami schen 
Computerprogramme erreicht, die auf analytischen Lósungen basieren. 
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Auch die in den Codes SING! und SING-S benutzten Randintegral- | [T] A. Ludwig, R. Krieg: An A ሚሚሚ E 
Verfahren bieten hier Vorteile. Fehler kónnen hier nur in den calculating eigenvibrations of thin circular cylindrical 


Randflüchen auftreten. Auch sehr feine Diskretisierungea, die shells. Journal of Sound and Vibration 74 (1981), im Druck 


bei Verwendung besonders effektiver Lósungsverfahren im Pro- 


— AÐA ሰቦ /^87 B. Göller: Dynamic Deformacions of thin spherical shells | 
based on analytical solutions. Journal of Sound and 


Vibration (1980), im Druck 
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Be Einleitung 





Die Qualifizierung von Rechenprogrammen für die Beschreibung des Brennstab- 
verhaltens beim Kühlmittelverluststórfall erfordert, neben der Durchführung 
von Einzeleffekt-Experimenten, auch Versuche mit nuklearem Brennstoff. Der- 
artige Versuche sind allerdings sehr aufwendig, so daß nur eine geringe An- 
zahl gefahren werden kann. Außerdem sind die Möglichkeiten der Versuchsin- 
strumentierung bei den in-pile-Tests (nicht allein wegen der Radioaktivität 
im nuklearen Core) stets sehr stark eingeschränkt. Daher muß zur Beurteilung 
der Aussagefähigkeit solcher in-pile-Tests eine anwendungsnahe Untersuchung 
der Eigenschaften der eingesetzten Versuchsinstrumentierung erfolgen, um Fehl- 


Interpretationen der nuklearen Versuche auszuschließen. 


Anlaß für die hier vorgelegte Untersuchung waren Diskrepanzen zwischen der 
Vorausberechnung der Hüllrohrterpersturen und den Meðergebnissen am LOFT- 
Reaktor ım Idaho National Engineering Laboratory, Idaho, USA ይፈ 

Die vorhergesagten Hüllrohrtemperaturen lagen, wie Abb. Í zeigt, über den 
ganzen Verlauf des simulierten Stórtalls auf einem Niveau von 2.5. etwa 

800 “e, wáhrend die gemessenen Werte durch ein nicht vorhergesagtes Benetzen 


der Hüllrohre durch das Kühlmittel meist um einige 100 K tiefer lager. Selbst 


die Maximaltemperaturen lagen bei der Messung noch um etwa 150 K tiefer als 


































bei der Rechnung. Es stellte sich daher die Frage, ob die vorhergesagten un- 
günstigeren Temperaturen auf einer Konservativitšt der verwendeten Computer- 
codes beruhen, oder ob bzw. zu welchem Anteil die bei LOFT eingesetzte Ther- 


moe lement-Instrumentierung der Brennstabhüllrohre für die festgestellte Dis- 


krepanz verantwort!ich sein kann / 3 /. 


Der Frage der Aussagefáhigkeit der Thermoelementanzeigen hinsichtlich der 
tats ichlich vorhandenen Hüllrohrtemperaturen unter den Stórfc' lbedingungen 
war selbstverständlich (vie es sich für gut vorbereitete Versuche versteht) 
bereits vor dor Durchführung der Experimente nachgegangen worden. Der expe- 
rimentelle Nachweis zu dieser Frage konnte aber zur Zeit der Ve rsuchs vorbe - 
reitung der nuklearen Tests nur sehr unvollkommen geführt werden, weil nir- 
gendwo geeignete Versuchseinrichtungen und Brennstabsimulatoren zur Verfü- 
gung standen. In dieser Frage eróffneten sich erst nach .'r Durchführung der 
nuklearen Tests neue MOglicukeiten, insbesondere auch mit der Ve rsuchs anlage 
COSIMA im Kernforschungszentrum Karlsruhe. 


2. Instrumentierung der Hiillrohre bei LOFT 








Bei der gekoprelten Untersuchung des Brennstabverhaltens unter Blowdown-Be- 
dingungen im Reaktor zur Messung cer Hiillrohrtemperaturen lassen sich nur 
Therroelemente verwenden, wobei diese nur auf die Brennstabhülle aufgesetzt 
werden kónnen. Um so wenig wie mglich mechanisch oder thermisch das Stabin- 
nere unc die Hille zu stören, wurde für LOFT die in Abb. 2 gezeigte Thermo- 
elementanordnung gewáhlt. An jedem instrumentier*en Brennstab sind 4 The rmo- 
elemente befestigt, welche Anzeigen für 4 verschiedene axiale Positionen lie- 
fern. Aus Symmetriegründen sind im beheizten Stabteil mit Hiife von Dummies 


alle 4 MeBleitungen über die ganze Stablánge durchgezogen. 


Die Anordnung der so instrumentierton 5rennstábe über den Reaktorkernquer- 
schnitt ist in Abb. 3 dargestellt. Als typisch für die LOFT-Instrumentierung 
erweist sich die Anordnung dieser Brenastübe in Gruppen zu je 5 Stáben. In 
dieser dichten Anordnung besteht bereits die Móglichkeit einer Rückwirkung 
der Instrumentierung (über die veränderten Wármeübertragungseigenschaften 
der Brennstábe mit den rippenartigen Ti ^rmoelemert-Mantelrohren) auf die 


sich in den Unterkanálen lokal einstellenden Kühlmittelzustände. 


In Abb. 4 und 5 sind Meßergebnisse für die LOFT-Transientez 12-2 und L2-3 


gezeigt, die mit dieser The rmoe lement-Inst rumentierung gewonnen wurden. In 


Abb. 4 erkennt man bei etwa 6 Sekunden ein scharfes Quenchen (Wiederbenet- 
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zen), das nach Abb. 5 über die ganze Core-Hóhe in Erscheinung tritt. Weiter- 
hin zeigt Abb. 5, daß im unteren und mittleren Kembereich die Brennstabtem- 
peracuren über längere Zeit wieder auf hohe Werte kommen, während im oberen 
Kernbereich noch dreimal ein Dryout mit anschlieSeadem Wiederbenetzen auf- 

tritt. Diese unerwartet stark zur Wirkung kommenden Wiederbenretzungserschei- 
nungen sind als wesentliche Ursache für die erwühnr. große Abweichung zwi- 

schen Vorausberechnung und Experiment anzunehmen. Es stellt sich die Frage, 

oh es sich bei dieser Wiederbenetzung überwiegend um einen Effekt handelt, 
der für die bei LOFT angewandte The rmoe lement-Inst rumentierung typisch ist, 
oder ob diese Wiederbenetzung auch bei den nicht instrumentierten Brennstá- 


ben auftreten sollte und von den konservativ (= pessimistisch) rechnenden 


Computercodes nur nicht wiedergegeben wird. 


3. Experimente mit COSIMA 





Zur Untersuchung dr Frage, ob die Wiederbenetzungserscheinung ein fiir die 
LOFT-Iostrumentierung typischer Effekt ist, bot sı.a die Versuchsanlage 
COSIMA im KfK an / 1 7. COSIMA ist für die Aufgabe konzipiert, out-of-pile 
das Wechselspiel zwischen mechanisch-thermischem Brennstabverhalten einer- 
selts und den thermohydraulischen Randbedingungen beiz Blowdown andererseits 
zu untersuchen mit dem Ziel, die Entwicklung von Computerprogrammen zu un- 
terstützen. Zu diesem Zweck wurde für COSIMA ein elektrisch beheizter Brenn- 
stabsimulator entwickelt, dessen Eigenschaften den nuklearen Brennstáben 

sehr nahe kommen (Original-Hüllrohr, Stableistung, Wármespeicherung, Gasspalt, 
Gasinnendruck, Pellet-Relocation). Weil hierbei das mechanische Verhalten der 
Hülle ein dominierender Effekt ist, mußte zur Messung der Hüllrohrtemperatur 
die berührungslose (und damit stórungsfreie) optische Messung mittels Pyro- 


metern herangezogen werden. 


Da es sich in COSIMA ohne gro&en Aufwand ermöglichen ließ, auch LOFT-typisch 
instrumeutierte COSIMA-Brennstabsimulatoren einzusetzen (Typ SIM), bot diese 
Versuchsanlage nahezu ideale Bedingungen zur Untersuchung des gestellten 
Problems. Das Interesse auf der amerikanischen Seite führte dazu, daß vor 
ungefáhr 1 1/2 Jahren zwei CJSIMA-Brennstabsimulatoren bei EG&G in Idaho mit 
Unterstützung durch die NRC mit einer LOFT-typischen The rmoe lement- Instrumen- 
tıerung versehen wurden. Inzwischen wurden mit diesen Simulatoren und 10 
nicht instrumeutierten Stáben insgesamt über 50 Blowdown-Versuche in COS IMA 
im Zusammenhang mit dem Nachweis des Thermoelementeffek tes durchgefiihrt. Ein 


Teil dieser Versuche diente dabei der ¡terati ven Einstellung geeigneter Steuer- 















Programme für den Blowdown-Ablauf und der Untersuchung der für die Reprodu- 


zierbarkeit der Messungen wesentlichea Stabparameter. 


3. |] COSIMA-Anlagenbeschreibung 





Über die COSIMA-Versuchsanlage war bereits auf dem PNS-Jahreskolloquium 1978 


(vgl. KfK 2770) ausführlicher berichtet worden / 17, so daß hier nur eine 


gekürzte Fassung der Anlagenbeschreibung gegeben vird. 


In Abb. 6 ist ein Schema der Versuchsanlage COSIMA (* conroled single rod 


blowdown simulation experiment Karlsruhe) dargest lit. Bei den Versuchen be- 
steht zunáchst die Aufgabe, cie stationären 
nachzubilden, 


Bedingungen des Reakcorbetriebes 
von denen aus der Blowdown gestartet wird. Daher bestebt die 

Anlage in ähnlicher Weise wie ein Druckwasserreaktor aus einem Wasserkreis- 
lauf mit Druckhalter, der den Systemdruck bestimmt. In einem elektrisch be- 
heizten Hei ßwasserkessel befindet sich auf 


der U 


geheiztes Wasser, welches mittels 
uwálzpumpe bis zum Erreichen stationárer Temperaturbedingungen im Kreis- 
lauf umgepumpt wird. 


Etwa 30 s vor der Blowdown- Aus lösung werden die kandheizstäbe und der zentral 


angeordnete Versuchsbrennst ab eingeschaltet. Durch die Wärmezufuhr zum Kreis- 


laufwasser bildet sich innerhalo der Teststrecke ein von unten nach oben an- 


steigendes Temperaturprofil aus. Kurz bevor das in den Hei Bwasse rkesse! 


fließende heißere Wasser dessen unteres Ende erreicht, muß die Blowdown-Aus lö- 


sung erfolgen, damit die gewollte Tlemperat:-schichtung in der Teststrecke er- 


halten bleibt. 


Der Betriebsdruck des Kreislaufes Leträgt etwa 160 38ኛ, die Betriebstempara- 


tur ist zum Zeitpunkt der Blewdown-Auslösung im oberen Kühlmittelplenum etwa 


A " er - - 
320 C. Die Anlage erstreckt sich über eine Hóhe von 18 m (Abb. 7). 


Die Teststrecke (vgl. Abb. 6) aat eine beheizte Länge von nur 1,5 m, weil an 


5 


den elektrisch beheizten Brennstaosimulator im Zentrum der Teststrecke beson- 


dere Anforderungen gestellt werden, die eine verkürzte Ausführung erforder- 


lich machen. Zum einen soll der Brennstab eine "bláh fahi ge” Zircaloy-Hülle 


erhalten und mit Innendruck beaufschlagt sein, und zum ande ren soll der Brenn- 


stoff hinsichtlich der im Reaktorstab gespeicherten Wirme, der Brennstoff- 


geometrie und schließlich auch der auf die Stablánge bezogenen Heizleistung 


nachgebildet werden. Diese Bediagungen werden mit dem in Abb. 8 gezeigten 


Stabkonzept erfüllt. Die maxima] möglıche Heizleistung beträgt im hoch be- 





zurück- 
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3.2 Beschreibi 





Für die Durch! 
verschiedene ] 
in Anlehnung a 
te ausgebildet 
Druckwasserre a 
sentlichen nur 
moe lement-Unte 
¡le nach etwa 

mildert wurd 
d Unterschie: 


erzielen, als 4 


Im Abb. 9 ist « 
COS IMA-Transier 
inst rumntierte 
dieser Vergleic 


geworfen hatte, 





leren Abschnitt des Brennstabsimulacors von 0,5 m Linge ca. 


cm. In Abb. 8 ist auch ein Pyrometer zur berührungslosen Messung 
temperatur dargestellt. 


"Auslösung wird durch das rasche Schließen der Schnellschlußven- 
1, wodurch die Teststrecke vom Kreislauf getreant wird. Mit einer 
tógerung werden sodann die Schnellóffnungsventile der beiden Aus- 
gen geöffnet, so daß der in der Teststrecke unter Druck stehende 
ı beheizte Wa»serinhalt in deu Kondensationsbehálter ausstrómt. 
men wird mittels der beiden in den Ausblaseleitungen befindli- 
¡teueventile in vorherbestimmter Weise gezielt beeinflußt. 


wird auch die Beheizung der Wandheizstäbe ausgeschaltet und die des 


lators zurtickgesteuert. 


Vorgang ist beendet, wenn der Druck ¿n der Teststrecke auf den 
s Meaktorsicherheitsbei.ilters abgesunken ist. In COSIMA ist es 
nacn dem Schließen der Auslaßventile durch Öffnen einer Ver- 
eißwasserkessel die Teststrecke mit heißem Wasser rasch zu flu- 
all ist dann von Interesse, wenn der gegen Ende des Blowdown er- 
nd des Simulators "eingefroren" werden muß, falls bei vorhande- 
L die Hüllrohrtemperatur noch hoch ist. 


ms der Blowdown-Transienten 





ührung der Untersuchungen zum Termoei2menteinflu8 wurden zwei 
"ansienten-Typen verwendet. Die áltere dieser Transienten wurde 
a eine zuvor in COSIMA sehr háufig gefahrene Blowdown-Transien- 
. Hiermit wurden die Folgen eines cold-leg 2F-Bruches bei einem 
ktor vom Typ Biblis A simulier:. Die Transiente umfaßt im we- 
der ersten Temperaturpeak zu Beginn des Blowdown. Für die Ther- 
rsuchungen wurde diese Transiente dahingehend abgewandelt, daß 

7 s nach Blowdown-Auslösung wieder einsetzende Kernkühlung etwas 
5. Das Ziel dieser Änderungen war, eine stärkere Differenzierung 
de im Quenchverhaiten der Stäbe mit und -ane Thermelemente zu 
lies mit der ursprünglichen Transiente möglich gewesen wire. 


las verwendete Steuerprogramm der Blowdown-Ventile für die sog. 

ite (= COSTRA) gezeigt, die zunächst für den Vergleich zwischen 

m und nicht instrumentiertem Stab eingesetzt wurde. Nachdem 

h einige neue Fragen zur Reproduzierbarkeit der Messungen auf- 
wurden mit COSTRA auch “ufangreiche Versuchsserien zum Stab- 
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verhalten bei aufeinanderfolgenden Blowdows mit demselben Stab und zur Repro- 
duzierbarkeit der Messungen mit Stäben aus unterschiedlichen bzw. gleichen 
Herstellungsserien gefahren (vgl. náchstes Kapitel). 


Das Steuerprogramm in Abb. 9 ist daraufhin ausgelegt, sofort zu Beginn des 
Blowdown-Ablaufes in der Teststrecke eine Strómungsumkehr zu bewirken, der 
nach einer kurzen Phase eines schwachen positiven Kerndurchsatzes ein lánge- 
rer Abschnitt der Strómung:stagnation folgt. In dieser Phase bildet sich im 
Kernbereich von oben her eine "Blase" trockenen Dampfes aus (dry out), die 
durch die nach etwa 7 s wieder einsetzende negative Kerndurchst römung (we- 
gen des Öffnens des unteren Abblaseventils) nach unten ausgewaschen wird. 

In dieser Phase erfolgt ein mehr oder weniger heftiges Quenchen des Brenn- 
stabes. 


Nach dem Vorliegen hinreichender Informationen über die LOFT-Transienten 
L2-2 und L2-3, es handelte sich im wesentlichen um den mit RELAP errechneten 
Kernmassendurchsatz in Abhángigkeit von der Zeit und um den Druckverlauf im 
Kern, konnte fiir COSIMA ein zweites Stewerprogramm iterativ erstellt werden, 
mit dem die Verháltnisse dieser LOFT-Transienten angenáhert nachgebildet 
wurden. Nach einer ersten Vorausberechnung des Steuerprogramms wurde die 


weitere Iteration experimentell-empirisch vorgenommen. 


Als ausreichen: Annäherung an die LOFT-Transiente L2-3 ergab sich das in 
Abb. 10 gezeigte Steuerprogramm für COSIMA; diese Transiente wird im fol- 
genden mit LOFTRA bezeicht. Trotz wichtiger Unterschiede zwischen der An- 
lagen LOFT und COSIMA besteht eine hinreichende Áhnlichkeit zwischen der 
Vorlage L2-3 und LOFTRA, um Aussagen ‘iber die miglichen Auswirkungen der 


LOFT-The rmoe lementinst rumentie rung bei den LOFT-Transienten zu machen. 


Ein wesentlicher Unterschied zwischen LOFTRA und COSTRA besteht darin, daß 
LOFTRA etwa die zwei fache Zeitdauer in Auspruch nimmt und dementsprechend die 
Abblaseventi le weniger weit geóffnet werden. Ein weiterer wesentlicher Un- 
terschied besteht im Verlauf des angestrebten Kernmassenstromes; nach der 
ersten Strómungsumkehr (typisch für einen großen cold-ieg-Bruch, noch ähn- 
lich wie bei COSTRA) Bibt es bei LOFTRA anstelle nur einer jetzt drei Stag- 
nationsperioden. Dazvischen wird zunáchst eine Aufwártsst römun g eingestellt 
(bei COSTRA war es eine Abwärtsströmung), darauf folgen nach Stagnation eine 
Abwárts- und nach emeuter Stagnation wieder eine Aufwärtsströmung (diese 


Vorgänge sind bei COSTRA nicht vorhanden). 


3.3 Reproduzierbarkeit der Mecsungen 





Die Ergebnisse der Messungen mit instrumentierten und nicht instrumentierten 
Stäben standen zunáchst im Widerspruch zur Erwartung. Dies betraf auch die 
Vergleichbarkeit der Messungen mit verschiedenen Simulatoren. Um den Vergleich 
verschiedener Stábe zu umgehen, wurde nach dem dritten Blowdown mit dem ersten 
instrumentierten Stab von diesem die Instrumentierung entfernt, so daß wieder 
ein nicht instrumentierter Stab entstand. Die damit in den folgenden Versuchen 
erreichten Hüllrohrtemperaturen zeigten jedoch anstelle der vorhergesagten Er- 
hóhung eine drastische Absenkung. Es handelte sich offenbar um eine Fehlein- 
schátzung des Brennstabverhaltens. 


Die Ursache dieser Fehleinschätzung wurde in den folgenden umfangreichen Me&8- 
reihen aufgeklárt. Sie bestand in der bisher ungenügenden Kenntnis des Stabver- 
haltens mit verformbarem Hüllrohr, dem inneren Gasspalt, sowie den beim Auf- 


he:zen zerspringenden keramischen Pellets im Simulator. 


bci den nun folgenden Messungen wurden diese verschiedenen Beeinflussungen sehr 
genau unter Kontrolle gehalten. Um die Vergleichbarkeit verschiedener Simulato- 
ren untereinander zu sichern, wurden nur noch Graphitheizstábe, Keramikpellets 
und Hüllrohre aus denselben Lieferungen (Chargen) verwendet, wie sie für die 
beiden instrumentierten Stähe zum Einsatz gekommen waren. Weiterhin wurde der 
Einfluß der Hüllrohrdeformation dadurch ausgeschaltet, da8 durch eine entspre- 
chende Wahl der Stableistung und des Stabinnendruckes die sich siss iienden 


Kriechverformungen erer etwa 10 um gehalten wurden. Dies wachte auch das Ein- 


führen des bereits im vorigen Abscanitt erwáhnten Zwangsquenchens der Brennstá- 
be am Ende des gesteuerten Blowdown-Ablaufes erforderlich, da ohne dasselbe 

in Teilbereichen der Stábe nach Blowdown-Ende unerwünschte und (wegen der un- 
kontrollierten Temperaturverláufe) nichtreyroduzierbare Eüllrohrverformungen 


auftraten. 


Bei der Durchführung von Blowdown-Serien gleichartiger Blowdowns mit jeweils 
nur einem Brennstab wurde gefunden, daß sich die Eigenschaften des Graphit- 
heizers so gut vie nicht veründerten, wührend die erreichten Hüllrohrtempera- 
-uren von Versuch zu Versuch systematisch und (mit einer gewissen Streuung 
der Daten) reproduzierbar kleiner wurden. Dieser Effekt ist beispielhaft in 
Abb. 11 dargestellt. Er erklärt sich durch eine "recolation" der von Ve rsuch 
zu Versuch in zunehmendem Maße zerspringenden Keramikpellets. Diese Zunahme 








ሀ 
— 


der Anzahl der Bruchstücke der Ringpellets konnte anhand von Stichproben an 


verschiedenen Brennstabsimulatoren einwandfrei nachgewiesen werden. 


Die Wirkung dieser Pellet-Relocation ist vor allem in einer Abnahme des wirk- 
samen Gasspaltes zwischen Pellet und Hüllrohr zu sehen. Hierdurch werden bei 
gleicher Heizleistung des Graphitheizers in der Stationürphase vor der Aus- 
lósung der Blowdown-Transiente die Temperaturen im Heizstab und im Pellet im 
mer mehr verringert, so da8 demeutsprechend auch die im Stab gespeicherte 
Wärmemenge beim Eintritt in den Blowdown w;rringert ist. Dies wirkt sich auf 
die Hüllrohrtemperaturen beim Blowdown schr wesentlica aus.  - Aus dieser Er- 
gebnis kann auch ersehen werden, daß die Fähigkeit der Wärmespeicherung ähnlich 
dem Nuklearstab für die out-of-pile-Simulation von Blowdown-Transienten von 
entscheidender Bedeutung ist. 


Die Versuche hatten nun klar gemacht, daß zum Vergleich der mit unterschiedli- 
chen Brennstabsimulatoren gewonnenen Ergebnisse nur solche herangezogen werden 
kónnen, bei denen jeweils der neue Simulator möglichst identisch in seinen Ei- 
genschaften ist mit den Vergleichsstäben, und bei denen die einzelnen Stábe 
eiue miglichst gleichartige Geschichte durchlaufen haben. Das bedeutet, daß 
bereits im Laufe des ersten Blowdown ein Stab mit Thermelement-Inst rumentie- 
rung (wegen der anderen Hiillrohrtemperaturen) sich ein wenig anders entwickelt, 
als ein Stab ohne Thermoelemente. Daher sind streng genommen jeweils nur die 
ersten Blowdown-Ergebnisse für jeden :tab vol! vergleichbar. Es ist allerdings 
nützlich, auch die Ergebnisse weiterer Blowdowns zu Kontrollzwecken heranzu- 


ziehen. 


Hinsichtlich der Reproduzierbarkeit der Hüllrohrtemperaturen verschiedener 
Simulatoren beim jeweils ersten Blowdown unter der Voraussetzung "identischer 
Stabfertigung" ist in Abb. 12 ein Beispiel gezeigt. Die Streuung der Versuchs- 
werte liegt unter etwa : JO K. Hierbei ist die Meßgenauigkeit der Pyrometer 
(*10 und -30 K) bereits enthalten. 


Trotzdem sind die in Abb. 13 gezeigten Kurven nicht ganz als typisch für die 
Reproduzierbarkeit der COSIMA-Me&ergebnisse zu bezeichnen, sondern zunächst 
ne; für die Reproduzierbarkeit der COSIMA-Versuche als Ganzes. Bei der kompli- 
zierten Steuerung der Versuche und bei der Leistungszufuhr zum Brennstabsimu- 
lator treten námlich gelegentlich kleinere hardwarebedingte Abweichungen auf, 
die erst nach der Versuchsdurchführung erkannt und quantitativ festgelegt wer- 
den kónnen. Wenn die Steuerkurven von verschiedenen COSIMA-Experimenten tat- 


sáchlich übereinstimmen, dann ergibt sich auch eine wesentlich kleinere Streu- 








ung der Hüllrohrtemperaturen als oben angegeben. Daher i 
barkeit der Meßdaten in COSIMA als sehr gut zu bezeichie 


Bei den im folgenden durchgeführten Vergleichen wurde au 
-erium der Übereinstimmung der tatsächlichen Steuerdaten 
chenden Versuchen besou«e rs geachtet. Daher kónnen an di 
16 gezeigten Ergebnisbeispiele deutlich engere Maßstäbe 
es die Reproduktionsversuche nach Abb. 12 nahe legen. 


4. Versuchsergebnisse 


Im folgenden werden die Versuchsergebnisse fiir die Trans: 
LOFTRA getrennt dargestellt, da sich die beiden Transient 
3.2 diskutiert, in ihren Eigenschaften sehr weitgehend ur 


Alle Versuchsergebnisse wurden einer weitgreifenden Analy 
möglichen Stórung.n und Unsicherheiten der Aussage unterz 
ein umfangreiches Datenmaterial eingehend untersucht und 

Diese Vergleiche betreffen im wesentlichen die Ventilstey 
nach den Versuchen tatsächlich vorliegen, die Drucke und 

Teststrecke, die Hüllrohrtemperaturen und die Stabinnendr 
stromverláufe in den Abblaseleitungen, die Wasserinhalce 

Plena, die Strukturtemperaturen und die Heizleistungen. D 
Hinweise auf das Geschehen der Teststrecke gewinnen, die 
meSbare Größen betreffen (wie z.B. Strómungszustand in de 
diese Weise ist eine recht weitgehende Versuchsinterpreta 
eine tragfáhige Basis für die Aussagen über den Einfluß d 
Instrumentierung auf die gemessenen Temperaturen gibt. 


Wie bereits erwähnt, wurden für unsere Untersuchungen 2 Bi 
vom Typ SIM mit der LOFT-Thermoelement-Instrumentierung ve 
te mit der TE-Instrumentierung zu den Transienten COSTRA ኒ 
eine Versuchsserie gefahren werden. Weil Sich die volle Ve 
Versuchen mit verschiedenen Stáben nur fiir den jeweils ere 
dem neuen Sizulator) ergibt, sollen im folgenden nur Beisp 
veils ersten Stab-Blowdown gezeigt werden. Auf diese Weise 
fangreichen Versuchsmaterial zu jeder Transientenart nur z 
zuziehen. 


In der COSIMA-Testst recke sind, verteilt über die beheizte 
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Stablánge von 





500 mm, (schraubenlinienförmig) acht Pyrometer eingebaut, die von unten nach 
oben von | bis 8 durchnumeriert sind. In Stabmitte, zwischen den Pyromete rme 8- 
stellen IN 4 und IN 5 ist ein DWR-Abstandshalter eingebaut, der beim Vorlie- 
gen von Nebelstrómung einen merklichen Einfluß auf die lokalen Kühlungsbedin- 
gungen stromabwärts ausübt. Etwas oberhalb und unterhalb der oeheizten Stab- 
lange sind ebenfalls derartige Abstandshalter angeordnet, die bei der entspre- 
chenden Strémungsrichtumg neben dem Einfluß des Leistungsprofils ebenfalls 
einen Einfluß auf die Würmeübergangsbedingungen an den benachbarten Pyrometer- 
zeßstellen IN Í und IN 8 ausüben. 


Die LOFT-typische TE-Instrumentierung umfaßt A Thermoelemente, die bei den 
COSIMA-Stäben axial genau auf die Höhe der geradzahligen Pyrometermeßstellen 

IN 2, 4, 6 und 8 positioniert wurden. Dagegen sind die Pyrometermeßstellen 
gegenüber den auf gleicher Höhe sitzenden Thermoelementen azimatal nur 45° 
versetzt angeordnet, so daß die Pyrometer, optisch wnbeeinflußt durch die Ther- 
moelemente, knapp 4,3 mm neben diesen die Hüllrohroberflächentemperatur erfas- 


Für die folgenden Vergleiche typischer Temperaturverláufe bei beiden Tran- 
sientenarten wurden jeweils die Pyrometermeßstelleu IN 4 (direkt unterhalb 
dem mittleren Abstandshalter) und IN 6 (etwa Mitte obere Stabhálfte), mit den 
zugeordneten Thermoelementen T 22.2 und T 23.1 herangezogen. Die Meßkurven fiir 


den instrumentierten und den nicht instrumentierten Stab werden jeweils zusammen- 


geplottet, um die Ergebnisse unmittelbar und anschaulich vergleichen zu können. 


4.1 Ergebnisse mit der Transienten COSTRA 





Die Abbildungen i2 und l4 zeigen die ausgewählten Beispiele aus den COSTRA- 
Vergleichsmessungen mit (Versuch V93) und ohne (Versuch V108) Thermoeiemente 
am Brennstabsimulator. Aus den zur Verfügung stehenden Versuchen ohne 
TE-Instrumentierung wurde hierbei derjenige herausgesucht, der bzgl. der 
Steuerdaten und Heizleistungsverläufe praktisch identisch mit V93 ist. Auf 


diese Weise sind die gezeigten Kurwen ohne Korrekturen vergleichbar. 


Es muß daran erinnert werden, daß die Pyrometer erst oberhalb einer Anzeige 
von knapp 810 K auswertbare Meßsignale lieferr, während die Thermoelemte na- 
türlich den Bereich der Fluidtemperatur noch voll erfassen. 


Folgende Aussagen wurden mit Hilfe der erwähnten Trendanalysen erhártet und 
werden durch die Abb. 13 und i4 für die Transiente COSTRA illustriert: 





during Steam Condensation. Nucl. Engrg. Des.(1981) im Druck 


- Für den nicht instrumentierten Stab liegen die Hüllrohrtemperaturen im 


allgemeinen über denjenigen für den iustrumentierten Stab (Ausnahme kurz- 


zeitig im Bereich um 35 s). 


- Die Thermoelementanzeige liegt im allgemeiaen tiefer als die Pyrometer- 


anzeige für denselben Versuch. 


- Die kritische Heizflichenbelastung (Temperaturanstieg) wird durch die 


Thermoelemente verspätet angezeigt, wobei die Verzögerungen stark schwan- 


ken und bis über | s betragen. 


- Die Peak-Cladding-Temperaturen werden von den Thermoelementen bis um 


ca. 30 K zu niedrig gemessen. 


- Bei kleinem Wirmeiibergang (Kühlung mit trockenem Dampf) verschwinden die 


Unterschiede zwischen Thermoelement- und Pyrometeranzeige. 


- Das Wiederbenetzen der Hüllrohre wird durch die TE-Instrumentierung 


stark begünstigt, wobei im selben Versuch die Thermoelemente vor den 


Pyrometern das Quenchen anzeigen. 


4.2 Ergebnisse mit der Transienten LOFTRA 





Die Abbildungen 15 und 16 zeigen die ausgewählten Beispiele aus den LOFTRA- 
Vergleichsmessungen mit (Versuch 126) und ohne Thermoelemente (Versuch 129) 
am Brennstabsimulator. Es handelt sich um dieselben MeCpositionen, wie bei 

den beiden Abbildungen zuvor. Als Vergleichsversuch ohne TE-Instrumentierung 
stand hier nur ein Versuch (V129) zur Verfügung, der allerdings hinsichtlich 


des Versuchsablaufs (Stewerdaten, Heizleistung) mit Versuch 126 gut überein- 


stimmt. 


Die durchgeführten Trendanalysen führen hier zu folgenden Aussagen, die im we- 


sentlichen durch die gezeigten Abbildungen 15 und 16 illustriert werden: 


- Zu Beginn der Transienten liegen die Hüllrohrtemperaturen (Pyrometeranzei- 
ge) für den nicht instrumentierten Stab höher als für den Stab mt Thermo- 
elementen. (Nicht im Bild zu sehen ist, daß die Temperaturunterschiede im 


ersten Temperaturpeak bei ca. 29 s lokal bis 150 K betragen.) 


- Mit den Thermelementen tritt die kritische Heizflächenbelastung verzö- 


gert auf. 





- Im Bereich mittlerer Wärmeübergänge, vermutlich im Zusammenhang mit Ne- 
belkühlung, treten besonders im ersten Yersuchsdrittel am instrumentier- 
ten Stab über viele Sekunden Dauer Thermoelementanzeigen auf, die um 100 


bis 150 K unterhalb der Anzeige des unmittelbar azimutal benachbart mes- 


senden Pyrometers liegen. 


- Der instrumentierte Stab neigt in der oberen Stabhálfte sehr stark zum 


Wiederbenetzen mit kurzzeitigen DNB-Erscheinungen. 


5. Diskussion des Thermoelementeinflusses 





Mit den beiden Transienten COSTRA und LOFTRA werden Bedingungen erreicht, 


die die beim Blowdown bei nicht oder nur wenig verformenden Hiillrohren mög- 


lichen Wärmeübergangsmodes zu einem wesentlichen Teil abdecken. Daran wird 


auch durch den Umstand nichts geändert, daß bei COSIMA die voll beheizte 


Stablánge nur 500 mm beträgt. 
Die im Experiment vorgefundenen Einflüsse der The rmoe lement- Inst rumentierung 


auf die jeweiligen Temperaturanzeigen bestátigen und quantifizieren die qua- 


litativen Erwartungen: 


Durch die LOFf-typische TE-Instrumertierung wird in der Hauptsache die Hü'lrohr- 


oberfläche vergrößert (Kühlrippen). Da die Thermoeiemente in diesen Rippen sit: 


zen, zeigen sie nicht die Temperatur des Hüllrohres, sondern diejenige der "Rip: 


" an. Weil jedoch diese "Rippen" nicht über ihre ganze Lange mit dem Hüllrohr 


pe 
verschweißt sind, neiren die thermisch schwach gekoppeiten Teile zur Bildung 
frühzeitiger lokaler Benetzungsfronten, sobald dazu Gelegenheit besteht. Daher 


sind die maximalen Meßfehler (bezogen auf die wahre Hüllrohrtemperatur in der 


näheren Umgebung) stets dann zu erwarten, wenn das kühlende Fluid in form von 


nassem Dampf (Nebelkühlung) vorliegt. Das Auftreten von vielen lokalen Bent- 


zungsfronten ist auch die Ursache für die Neigung der instrumentierten Hüll- 


rohre zum vorzeitigen Wiederbentzen des ganzen Brennst abes. 


Bei den übrigen Wármeübergangsmodes sind die Abweichungen der verschiedenen 
gemäß viel kleiner und hier auch durchaus (unter 


Temperaturanzeigen erwartungs 
Bei der kritischen 


Berücksichtigung der schwierigen MeBaufgabe) tolerierbar. 


Wármebelastung macht sich die um etwa 50 Z vergrößerte Staboberfläche in einer 
g des Eintritts des DNB bemerkbar. In dieser Zeitspanne 
als es bei einem Brennstab ohne Thermo- 


kleinere reak-Cladding- 


deutlichen Verzögerun 
verliert der Brennstoff mehr Wärme, 


elemente der Fall wäre; dementsprechend werden auch 








Temperaturen erreicht. 


Werden zeitweilig kleine Wärmeübergangszahlen bei reiner Dampfkühlung erreicht, 
dann findet zwischen Hüllrohr und Thermoelementen ein Temperatursusgleich statt, 
so daß die Anzeigeunterschiede verschwinden. Gelegentlich kann der Dampf jedoch 
gegenüber dem Hüllrohr auch überhitzt sein. In diesem Fall werden in den Expe- 
rimenten gegenüber den Pyrometern erhöhte Thermoelementanzeigen vorge funden 
(die gezeigten Abbildmgen enthalten dafür kein Beispiel). 


Diese qualitative Diskussion zeigt, daß der vermutete Thermoelemeuteinflu8 im 
Einklang mit den Erwartungen steht. Eine quantitative Diskussion kann nur an- 
hand der vorgelegten experimentellen Daten durchgeführt werden, da bis heute 

keine Rechenprogramme existieren, mit denen der Einflu8 der Thermoelemente ua- 


ter Blowdown-Bedingungen quantitativ erfaßbar ist. 


6. Sch lußfolge rungen 
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Beim nuklearen Brennstab und bei den COSIMA-Brenustabsimulatoren geht die i 
Stab erzeugte Würme von den Keramik-Pellets über einen durch "Fuel Relocaticn" 
beeinflußten Gasspalt, sowie über das dünnwandige Cladding an das Kühlmittel 
über. Vom Standpunkt der Würmeübertragung war zu erwarten, daß bei den Vorgán- 
gen beim Blowdown ein merklicher Einfluß der LOFT-typischen Thermoelement- 


wei Instrumentierung auf die als Hüilrohrtemperaturen gedeuteten Temperaturanzei- 
\- | gm der Thermelemente auftritt. 
Bei den in COSIMA durchge führten Experimenten konnte gezeigt werden, daß die- 
ser Einfluß beim Vorliegen von Mebelkühlung und hinsichtlich des Auftretens 
der Wiederbentzung von entscheidender Bedeutung für die Versuchsinterpreta- 
tion sein kann. 
E: ሐሐ 
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Abb.1: LOFT- Versuch L2-3, Vergleich gemessener und gerechneter Heiflstabtemperaturen 
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Abb. 8: COSIMA - Simulationsbrennstab SIM und Pyrometer 
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Lichtleiter 
Schnitt durch die Teststrecke mit Pyrometer 
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Abb.7: COSIMA - Gesamtansicht 
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Abb.6: Schema der COSIMA- Anlage 
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